function xOut=sum(xIn,is)
  if(nargin<2)
    is = 1;
  end
  xOut.v = sum(xOut.v,is);
  xOut.d = sum(xOut.d,is+1);
  xOut=class(xOut,'der1');
